Ultrafast Nonlinear Optical Response of Strongly Correlated Systems: Dynamics in 

the Quantum Hall Effect Regime 
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We present a theoretical formulation of the coherent ultrafast nonlinear optical response of a 
strongly correlated system and discuss an example where the Coulomb correlations dominate. We 
separate out the correlated contributions to the third-order nonlinear polarization, and identify 
non-Markovian dephasing effects coming from the non-instantaneous interactions and propagation 
in time of the collective excitations of the many-body system. We discuss the signatures, in the 
time and frequency dependence of the four-wave-mixing (FWM) spectrum, of the inter-Landau level 
magnetoplasmon (MP) excitations of the two-dimensional electron gas (2DEG) in a perpendicular 
magnetic field. We predict a resonant enhancement of the lowest Landau level (LL) FWM signal, a 
strong non-Markovian dephasing of the next LL magnetoexciton (X), a symmetric FWM temporal 
profile, and strong oscillations as function of time delay, of quantum kinetic origin. We show that 
the correlation effects can be controlled experimentally by tuning the central frequency of the optical 
excitation between the two lowest LLs. 
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I. INTRODUCTION 

The properties of systems far from equilibrium and, 
in particular, the role of many-body and collective 
effects on the femtosecond and the nanometer scale 
present relatw&Ujr-unexplored frontiers of condensed mat- 
ter physics E&ffiffl Such problems are particularly chal- 
lenging in semiconductors, where the time intervals of 
interest are often shorter than the interaction time; 
oscillation periods of the elementary excitations!! 
Examples of well-established pictures for the interaction 
processes that need be revised in this regime include the 
semiclassical Boltzmann picture of point-like particles ex- 
periencing instantaneous collisions and-the thermal bath 
pictures of relaxation and dephasing. ooQ Even the no- 
tion of weakly interacting "quasiparticles" , a cornerstone 
of condensed matter physics, must be revisited when de- 
scribing the ultrafast nonlinear optical response^. 

Wave-mixing experiments are ideally suited for ex- 
ploring quantum coherence and collective and.-correla- 
tion effects in semiconductor nanostructuresBBa Time- 
dependent interactions and correlations dominate the 
FWM signal during negative time delays, where the Pauli 
blocking effects vanish.Era The treatment of such inter- 
actions within|-.the time-dependent Hartree-Fock (HF) 
approximation^ predicts an asymmetric temporal profile 
of the FWM signal.Eratl The negative time delay signal 
generated by mean field exciton-exciton interactions de- 
cays twice as fast as the positive time delay signal. The 
observation of strong deviations from this asymmetric 
HF temporal profile in undoped semiconductors was at- 
tributed to exciton-exciton correlations aa 

The importance of many-body effects in determining 
the time and frequency profile of the ultrafast nonlin- 



ear optical spectra may be traced microscopically to the 
coupling, via the interactions, of the one-particle density 
matrix that describes the optical polarization measured 
in the experiment to many-particki-aprrelation functions 
(e.g. the higher density matrices). Bcm The latter are fa< 
torized within the time-dependent HF approximation 
The correlation-induced fluctuations, described by the 
deviations from the factorized form, generate a new 
FWM signal, which can display a distinct time and fre- 
quency dependence as compared to the mean field sig- 
nal. Such correlation effects are most pronounced during 
time scales shorter than the characteristic times associ- 
ated with the interaction processes.oEl 

To describe the above non-equilibrium many-body ef- 
fects, one must use a controlled truncation of the infinite 
hierarchy of coupled density matrix or Green function 
equations. In undoped semiconductors, this hierarchy 
truncates iLanfi-adopts an expansion in terms of the op- 
tical fields oBqEJOIIj This is the case since (a) in the 
ground state, the conduction band is empty and the va- 
lence band is full, and (b) the Coulomb-induced coupling 
of the conduction and valence bands via, e.g., Auger-like 
processes is negligible: in the absence of optical fields, the 
numbers of conduction band electrons and valence band 
holes are independently conserved. In undoped semicon- 
ductors, the lowest electronic excitations of the ground 
state electrons are the high energy interband e-h pairs, 
which can adjust almost instantaneously to the dynam- 
ics of the photoexcited carriers .Ej The photo-excited e- 
h pairs then behave as quasi-particles with mutual in- 
teractions, while the ground state can be considered as 
rigid. In this case, the many-body nature of the system 
only affects the different parameters associated, with the 
band structure and the dielectric screening, £-3 and the 
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only Coulomb correlations that require consideration are 
dynamically generated by the optical excitation.tj The 
almost unexplored dynamics of strongly correlated sys- 
tems, whose ground state electrons interact unadiabati- 
cally with the photo-excited e-h pairs, raises very funda- 
mental questions. 

A widely used theoretical approach for treating 
the above many-body effects in undoped semiconduc- 
tors is thr,, "dypamics-controlled truncation scheme" 
(DCTS).iyiiiEl In this theory, the response of the semi- 
conductor is expanded in terms of the number of created 
e-h pairs. Importantly, the Coulomb interactions that 
contribute to a specified order in the applied field only 
occur between such e-h pairs. This is the case since there 
is the correspondence between the number of e-h pairs 
and the sequence of photon absorption and emission, and 
there are no carriers in the ground state to interact with 
the photoexcited carriers. The latter condition is not 
met however in doped quantum wells, where a correlated 
2DEGLi£ present in the ground state, and the DCTS fails 
there.EiJ A new method that extends the DCTS principles 
to systems with a strongly correlated ground state is re- 
quired. In a series of works we applied a theory based on 
a canonical transformation and time-dependent coherent 
states to study the case where the interactions between 
the photoexcited e-h pairs and the electron Fermi sea 
(FS) excitations dominate the coherent nonlinear optical 
response. BEjO 

In FS systems, the direct exciton-exciton interactions, 
which dominate the nonlinear response in undoped semi- 
conductors, are screened, and the nonlinear response is 
determined by the FS excitations. For resonant pho- 
toexcitation, the optical dynamics is dominated by_ jit. 
elastic electron-electron (e-e) scattering processes. Eat£l 
At low temperatures, the dephasing times close to the 
Fermi edge increase by a-few picoseconds, in agreement 
with Fermi liquid theory.ELS For below-resonance excita- 
tion, however, the dissipation processes are suppressed 
and coherent effects dominate. A novel dynamics of 
the Fermi Edge Singularity is then observecOEll, due to 
many-body correlations, of the photoexcited holes with 
the FS excitations BOEa 

In the absence of long-lived excitations, a many- 
particle system, such as a FS, interacts with the pho- 
toexcited e-h pairs almost instantaneously, i.e. during 
time scales shorter than the pulse duration. The sys- 
tem then behaves to first approximation as a thermal 
bath, and its interactions with the photoexcited carriers 
can be treated within the dephasing and relaxation time 
approximations. This is not the case however if the du- 
ration of the interactions-is comparable to or longer than 
the measurement times.EH In the latter case, the semi- 
classical instantaneous collision picture breaks down, and 
quantum mechanical interference effects lead to noniix^i-t. 
nential decay and non-Markovian memory effects. BulU'Ell 
To study dephasing in the above quantum kinetic regime, 
one must account for the time pjYolution of the coupled 
photoexcited carrier-FS systemBtll 



The change in the energy spectrum caused by a per- 
pendicular magnetic field restricts the phase space avails 
able for e-e scattering in doped quantum wells (QW).ca 
For strong magnetic fields, the Coulomb correlations 
are enhanced due to the suppression of the kinetic 
energy.a In the quantum Hall effect (QHE) regime,E3Efl 
long lived .piUacti.ve excitations dominate the 2DEG 
spectrum.t3'Ejl£3'L3 Recently, the first experimental stud- 
ies of the role of such collective excitations in. the, ultra- 
fast nonlinear optical dynamics were reported.ooEj The 
presence of low energy excitations and the strongly corre- 
lated ground state raise formidable theoretical difficulties 
for describing the dephasing dynamics of the 2DEG. 

We are interested in developing a theoretical frame- 
work for describing the ultrafast dephasing and the non- 
linear optical response of strongly correlated systems. 
Examples of systems of interest include modulation- 
doped semiconductor QWs, where different strongly cor- 
related ground states are realized in the QHE regime, and 
the ferromagnetic semiconductors doped with magnetic 
impurities. In the first part of the paper (sections II-V) 
we describe the third-order nonlinear optical response of 
a many-electron two-band system without assuming a 
HF or other specific ground state. In the second part 
(sections VI- VII) we study the role of the inter-LL MP 
collective excitations in the transient FWM spectrum of 
the cold 2DEG. Here we concentrate on filling factors 
close to v = 1, where the spin-| ground state electrons 
lead to ferromagnetic properties (QHE ferromagnet), and 
the excitaticp. s»|ectrum is governed by strong Coulomb 
correlations. c3'E3'L3 We consider photoexcitation with tr+ 
circularly polarized light, in which case only spin-| elec- 
trons are excited and the_MP collective excitations play 
the most important role.EH Our results explain the most 
salient qualitative features of the transient FWM spec- 
trum observed in recent experiments Ej 

Our theory applies to a two-band system described 
by a Hamiltonian that independently conserves the num- 
ber of conduction band electrons and valence band holes, 
e.g. the GaAs/AlGaAs QWsIj We describe the coupling 
to the optical field within the dipole approximation, and 
neglect any stimulated emission. We consider zero tem- 
perature, which is adequate for describing correlations 
that require thermal energies smaller than the excitation 
and interaction energies of the system in order to be ob- 
served. The third-order polarization calculated here is 
expected to describe the nonlinear optical signal when 
the photoexcited carrier density is smaller than the den- 
sity of the ground state electrons, in which case the cold 
2DEG correlations prevail. 

The outline of the paper is as follows. In Section || we 
set up the general problem and discuss the nature of the 
states that contribute to the optical spectra. In Section 
III we study the time evolution of the system, and in- 



troduce a decomposition of the photoexcited many-body 
states that allows us to classify the different interaction 
contributions. In Section IV we use the above decom- 



position to derive the equation of motion for the third 
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order nonlinear polarization, Eg . (37). The decompo- 
sitions introduced in Section |l| allow us to distinguish 
the coherent and excitonic effects from the incoherent 
effects, and separate out the factorizable from the cor- 
related nonlinear polarization contributions even in the 
case of a strongly correlated ground state, In Section |v| 
we discuss an example of a basis of strongly correlated 
states that can be used to obtain equations of motion for 
the correlation functions that describe the many-body 
effects. In Section.- jiO.-we-.derive a generalized average 
polarization model, 0He3'E3 which we use to identify the 
signatures of the collective 2DEG excitations in the time- 
dependent FWM spectra. The ground state correlations 
determine the interaction parameters in the equations of 
motion. We derive in the appendices a number of re- 
lations among such interaction parameters that are im- 
posed by t he e lectron-hole symmetry of the ideal 2DEG. 
In Section VII we present numerical results that describe 



the correlation-induced ultrafast dynamics predicted by 
the above model. We identify a number of interesting 
features in the time-dependent FWM spectrum, which 
arise from the propagation in time of the inter-LL MPs 
and their non-instantaneous interactions with the pho- 
toexcited excitons. We end with the conclusions. 



with cyclotron energy SI" (see Appendix [aJj.B The opti- 
cal transition operator X* is expanded in terms of inter- 
band e-h pair creation operators Jtj that we refer to as 
the exciton (X) operators from now on: 



X 



f = E 



'NiX 



(3) 



In the case of the 2DEG in a magnetic field it is conve- 
nient to introduce the LLi magnetoexciton states = 
Xj\0), where |0) is the ground eigenstate of the many- 
body Hamiltonian H , with full valence band and the 
2DEG at rest. The eigenvalue equation H\0) = de- 
fines the ground state energy as the reference point. In 
the ideal system, 
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(4) 



where iV; = N(l 



with N = L 2 /2irl 2 being the LL 



degeneracy, I the magnetic length, L the system size, and 



(5) 



II. PROBLEM SETUP 

We are interested in developing a comprehensive ap- 
proach to the problem of the nonlinear optical response 
in the case of photoexcitation from the valence to the 
conduction band. Within the dipole approximation, the 
coupling to Ahe optical field can be described by the 
HamiltonianO (7j = 1) 



fltot(t) =H- n£{t)& - nE*{t)X. 



(1) 



In the above equation, H is the "bare" many-body 
Hamiltonian, which describes the bandstructure effects 
and the interactions, £(t) is the applied optical field, X 
is the optical transition operator, and /i is the interband 
transition matrix element. In the case of a semiconductor 
QW containing a 2DEG in a magnetic field, the Hamil- 
tonian H has the formtl 

ff = 2£ ( + nS(i + l/2)]4 ii e w 
+ J2 Wit + 1/2) hl k /h^ % + V ee + V hh + V eh , (2) 

where E g is the bandgap, and V ee , V e h, and Vhh are, re- 
spectively, the e-e, e-h, and h-h interactions (see Ap- 
pendix ^) . The magnetic field splits the conduction and 
valence bands into discrete electron (e) and hole Qi) LLs, 
e-LLi and /i-LLi, where i includes both the LL index 



and the spin. 



is the creation operator of the LLi 



conduction band electron, with cyclotron energy and 
hi i is the creation operator of the LLi valence band hole, 



gives the filling of LLi in the absence of optical excitation. 
Note that the exciton states \Xi) are strongly correlated: 
they are created by the operator xj acting on the ground 
eigenstate of the many-body Hamiltonian H, which de- 
scribes the correlated electron gas at rest. From Eq. (H) 
we obtain the commutation relation 



[X i ,X]} = S ij 1 



AN; 



where the operator 



Nui, 



(6) 



(7) 



with (0|AAi|0) = 0, describes the number of photoex- 
cited carriers in LLi. 

The optical spectra are determined by the polarization 
of the photo-excited system, 



P(t) = li^m) = m]T y/NiPi(t), 



(8) 



where Pi are the average values of the exciton operators 

Pi(t) = (V|*#>. (9) 

The state \ip(t)) evolves from the state of the system prior 
to the optical excitation according to the Schrodinger 
equation for the Hamiltonian H to t(t). For zero tempera- 
ture, this initial state is the lowest many-body eigenstate 
|0) and describes all correlations in the absence of optical 
fields. 
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As in the theoretical approaches of Refs.|ll|,y_2] there 
is a one to one correspondence between the photon ab- 
sorption/emission and the e-h pair creation/destruction. 
Since an electron gas may be present in the ground state, 
we classify the photoexcited states in terms of the num- 
ber of valence band holes, i.e. the number of missing 
valence band electrons as compared to the ground state 
|0). We thus decompose the optically-excited state \tjj) 
as 



l^o) + |V>i> + 1^2), 



(10) 



where \ip n ) is the collective n—h photoexcited-state. The 
above holes interact strongly with the 2DEG.O Note that 
the states with n > 3 da not contribute to the third-order 
nonlinear polarization.D 

Substituting Eq. (|l^) into the Schrodinger equation for 
the Hamiltonian H to t(t) we obtain up to third-order in 
the optical field that 

id t \ip )-H\ipo) = -ft£*X\ip 1 ), (11) 

- H\*h) = -v£Xi\1>o) - v£*X\i> 2 ), (12) 

idM - H\fh) = (13) 

with initial condition \ip n (~ 00)) — <5 n) o|0), where the 
Hamiltonian H includes the degrees of freedom that lead 
to the dephasing. The physics of the above equations 
is clearly displayed: |-0o) is coupled to \ipi) by the de- 
struction of one e-h pair, \ip±) is coupled to \ip2) by the 
destruction of one e-h pair and to \ipo) by the creation of 
one e-h pair, and \i>2) is coupled to \4>i) by the creation 
of one e-h pair. Fig. [h shows the optical transitions that 
determine the FWM signal up to third order in the op- 
tical field. It is worth noting that, by retaining in the 
expansion Eq. (|l0|) states with higher h numbers, one 
can extend Eqs. (gl]-[ij|) to treat higher order nonlinear 
processes. 

Even if we restrict ourselves to the electronic degrees 
of freedom, the Hilbert space of states that determine 
the ultrafast nonlinear response of a doped QW is com- 
plicated. Strictly speaking it contains all the states that 
can be generated through the coupling of e-h pairs pho- 
toexcited in any of the QW subbands with all the excita- 
tions of the 2DEG: plasmons, magnons, incoherent pairs, 
etc. For the purpose of developing an intuitive picture 
of the important physical processes it is useful to first 
discuss qualitatively the ensemble of states that are most 
relevant to the problem at hand. 

For the experimental conditions considered in the sec- 
ond part of the paper, the most important 2DEG excita- 
tions are the collective inter-LL MP modes, which arise 
from the coherent promotion of a LLO electron to a higher 
LLOEj iSjich MP eigenstates are well approximated by 
the formE3 

\^)=J^Pjf(q)e{ +qyJ e k>j ,\0), (14) 

kjj' 

where |0) is the strongly correlated ground state and the 
amplitudes pjj> (q) are related to the LLj' — » LLj contri- 
bution to the density operator. Note that, similar to the 
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FIG. 1: Photoexcitation of the intermediate (a) 0-h, \ipo), 
and (b) 2-h, \tp2), states via the nonlinear optical processes 
that contribute to the FWM spectrum. To third order in 
the optical fields, the coherent emission of a k s photon in the 
FWM direction k s — 2fc2 — k\ is determined by the excitation 
of two e-h pairs by the optical field A?2, and the deexcitation 
of one e-h pair by the optical field ki. Although in a coher- 
ent FWM experiment we must begin and end the nonlinear 
excitation process with the system in its ground state, the in- 
termediate 0-h state does not need to be the ground state |0), 
but can contain electron gas excitations. The above optical 
transitions are assisted by Coulomb interactions, which lead 
to the correlations discussed in Section [II. 



exciton states \Xi), the above MP states are strongly cor- 
related. For the magnetic fields of interest, in the ground 
state |0), only e-LLO is partially filled with the 2DEG 
at rest, while all the h-LL states are empty (full valence 
band) . Since we focus on photoexcitation of the LLO and 
LL1 optical transitions, the main contribution to the op- 
tical spectra comes from the resonant LLO — > LLlpMPs, 
whose energy is close to the LLO — > LL1 energy.E§c3 

It is useful to make the junction with two domains 
well studied in the recent literature: photoexcited un- 
doped QW, and 2DEG in the QHE regime. One can 
distinguish between the excitations of two subsystems: 
(i) the QW interband excitations (with the 2DEG at 
rest), which consist of le-h, 2e-h, ■ ■ ■ pairs created in the 
different QW LLs, and (ii) the 2DEG excitations (with 
unexcited QW and full valence band), i.e. the 1-MP, 
2-MP, ■ ■ • states, etc. The ensemble of states that deter- 
mine the third-order nonlinear optical spectra can then 
be thought as consisting oiie-h pairs, I < 2, and n 2DEG 
excitations. For photoexcitation of the LLO and LL1 ex- 
citon transitions, the inter-LL MP provides a resonant 
coupling of the two LLs, since its energy is comparable 
to the LLO — > LL1 excitation energy. In contrast, the LLi 
exciton states with i > 3, the states with n > 2 MPs, and 
the continuum of incoherent 2DEG pair excitations anal- 
ogous to the ones in an ordinary Fermi liquidEij primarily 
contribute to the optical spectra via non-resonant pro- 
cesses. 

One can draw an analogy between the X-MP effects 
of interest here and the X-phoncm|-| m t or act i on effects 
studied in undoped semiconductorsflu'HOEirEJ However, 
there are some important differences. In the undoped 
system, the electronic operators commute with the collec- 
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tive excitation (phonon) operators, and the ground state 
correlations can be neglected. One can then expand the 
state \ip) in terms of a basis consisting of products of 
phonon wavefunctions times e-h pair two-particle wave- 
functions. In contrast, a MP is an electronic excitation 
(see Eq. (fL4|)), and its creation operator may not com- 
mute with other electronic operators. Pauli exchange ef- 
fects must then be considered, while, unlike for phonons, 
MPs do not strictly obey Bose statistics. Importantly, 
one must treat the strong correlations of the ground state 
electrons. Issues such as the above complicate the use of 
a simple basis to calculate the nonlinear optical response 
of the 2DEG. In Section |v| we discuss an example of a 
strongly correlated basis set that can used to address the 
above issues. An important advantage of this particular 
basis is that it facilitates the development of a simple 
model that describes the most salient dynamical features 
of the ultrafast nonlinear optical spectra. 

III. TIME-DEPENDENT INTERACTION 
EFFECTS 

In this section we consider the time evolution of the 
coupled photoexcited carrier-2DEG system that leads to 
the dephasing of the e-h polarization. We are mainly 
interested in dephasing due to electronic degrees of free- 
dom, and thus the distinction between the photoexcited 
carriers and the "bath" excitations is less clear as com- 
pared e.g. to the case of a phonon bath. We address this 
issue by separating out the excitonic contribution directly 
excited by the optical field (2DEG at rest) from the con- 
tribution of the excited 2DEG configurations (denoted 
by 2DEG* from now on) that lead to the dephasing. For 
this we decompose the \-h photoexcited state as follows: 

|^i)=^^ L |X i ) + |^) ) (15) 

i 

where is the {l-/i/2DEG*} contribution defined by 
the condition (Xj|^i) = 0, and the exciton amplitude 

Pt = (X^) = (0|*#> (16) 

reduces to the linear polarization to first order in the op- 
tical field. To describe the two-photon nonlinear optical 
processes in Fig. ([!]), we must consider, in addition to the 
X-2DEG interactions, the X-X and X-|^i) interactions 
during the optical transitions. For this we first separate 
out the total interaction contribution to the 2-h and 0-h 
intermediate states, l?/^) and |^o) respectively, and then 
identify the particular contributions due to the interac- 
tions among the above 1— h excitations that lead to the 
correlation effects: 

ii' 



where the state \XiXit) — XjXj, |0) describes two non- 
interacting Xs, and 

|Vo) = (o|V> |o) + W l ) , 

i 

where we have separated out the ground state contri- 
bution from the 2DEG* contributions by requiring that 
(0|^ nt ) = Wo) =0. 

The above decompositions are analogous to the cu- 
mulants introduced within^-the DCTS for the case of 
undoped semiconductors.uO Such cumulants were ob- 
tained by subtracting the factorized contributions from 
the many-body correlation functions. Note however that 
the method presented below also holds in the case where 
strongly correlated carriers are present in the ground 
state, as in the 2DEG case, where the assumptions of 
the DCTS break down. The decompositions of the states 
IV'™*), IV^"*)' and l^i) also allow us to separate out, in 
the equations of motion Eqs. (|i"Tl), (|l2|), and (|l3|), the 
source terms proportional to the optical field from the 
source terms proportional to the polarizations Pf 1 , which 
lead to different time dependencies. The photo-excited 
states l^o), IV'i)) and |^) describe correlated contribu- 
tions, whose physical origin will be discussed below. 

We now derive the equations of motion of the above 
photoexcited states, which we will use in the next section 
to derive the nonlinear polarization equation of motion 
and separate out the factorizable contributions. It is eas- 
iest to start with the 1-h time-evolved state. Eq. jl5| ) 
splits this state into excitonic (2DEG at rest) and {1- 
/i/2DEG*} parts, Pf(t) and $1) respectively. The state 
originates from the X-2DEG scattering during the 
time evolution of the photo-excited X. To describe such 
interactions, we consider the action of the Hamiltonian 
H on the exciton states \Xi). By subtracting all the ex- 
citon contributions, the state if can be expressed in 
the form 

H\X t ) = ni\Xi) - £ Vi,i\X v ) + \Yi) , (19) 

where 

fli = (Xi\H\Xi) (20) 

is the Xi energy, 

V Vi = -(X v \H\X i ) = V& (21) 

describes the Coulomb-induced coupling of the different 
Xs, and \Y) = Y \0), where the operator 

Y = [X,H] - Q t X t + J2 Vn'X>, (22) 

describes the interactions between Xj, and the rest of the 
carriers present in the system. 
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FIG. 2: Photoexcitation of the exciton state, and then 

of the X+MP state |Yi) via resonant X-2DEG scattering. 



As one can see by using the above equations, the state 
\Yi) is orthogonal to all exciton states \Xj), (Yi\Xj) = 0, 
and is therefore the {l-/i/2DEG*} state into which Xi 
can scatter by interacting with the 2DEG. For the exper- 
imental conditions of particular interest here, the most 
important contribution to \Yi), Eq. (A17), comes from 
X+MP states. To see this, let us consider the possible 
final scattering states of the LL1 exciton X\. Its LL1 
electron can scatter to LLO by emitting a LLO — > LL1 
MP, a process shown in Fig. ||. Since the MP energy is 
close to the e-LLO — > e-LLl energy spacing, the above 
scattering process is almost resonant. It therefore pro- 
vides an efficient decay channel of the LL1 exciton to a 
{1-MP + 1-LLO-c + 1-LLl-h} four-particle excitation of 
the ground state |0). All other allowed scattering pro- 
cesses are nonresonant. The X\ hole can scatter to LLO 
by emitting a MP, which leads to a {1-MP + 1-LLl-e + 
1-LLO-h} four-particle excitation. The latter state how- 
ever has energy that is significantly higher, by an amount 
of the order of ~ S!£+f^ , from that of the initial X\ state. 
Note that, as shown in Appendix in the electron-hole 
symmetric limit the X electron or hole must change LL 
during the scattering process. In the case of X , the LLO 
electron can scatter to LL1 by emitting a MP, so that 
X -> {1-MP + 1-LLl-e + 1-LLO-h}, or the LL1 hole 
can scatter to LLO, in which case X$ — ► {1-MP + 1-LL0- 
e + 1-LLl-h}. \Yq) is thus a linear combination of the 
same final states as \Yi). However, in this case the en- 
ergy of all final states is significantly higher than that of 
the initial state \Xq). Therefore, the decay of the LLO 
exciton is suppressed as compared to that of the LL1 
(or higher) exciton. Note that the distinction between 
resonant and nonresonant interaction processes is most 
pronounced when the inter-LL excitation energy, of the 
order of the cyclotron energy, exceeds the characteristic 
2DEG Coulomb correlation energy, ~ e 2 /l. Expansions 
in terms of the ratio of the above two energies are kncpju, 
to capture most of the 2DEG correlation effects 



We now describe the time evolution, to first order in 
the optical field, of the 1-h photo-excited state \ipi). The 
equation of motion for the linear polarization Pf 1 can 
be derived by projecting the exciton state (Xi\ to the 
truncated Eq. (p^) and applying Eq. ([il)|): 



id t Pt = ihlf - V u ,pt + Pt- ^S{t)Nl /2 . (23) 



The correlation function 



P L 



(24) 



discussed in Section |v|, describes the dephasing of Pf 4 
and screening effects. 

Substituting the decomposition Eq. ([l5]) into the 
Schrodinger equation Eq. (|12}) , and using Eqs. (|23| ) 
and (flS|), we obtain the equation of motion of the {1- 
/i/2DEG*} contribution 



E 



[ptm-pt^)] 



(25) 



The operator P~Y? — P^Xj also appears in the equation 
of motion of the 2h state. Its first term describes the 
scattering of A", with the 2DEG, while its second term 
compensates for the dephasing of Pf 1 and ensures the 
orthogonality = 0. 

We can perform a similar analysis of the time-evolved 
2-h state by separating out in Eq. (17) the contribution 
of the non-interacting two-exciton states |XjXi/)£3'ED 
This contribution describes the time evolution of the two 
Xs photo-excited by the optical field in the absence of in- 
teractions. However, the two Xs interact with each other 
as well as with the 2DEG, as described by the equation 

H\XiXi,) = (fij + (U^XiXi,) - J2 Vji>\XiXj) 
- VjilXjXi.) + \XiY v ) + \XfYi) + \B U ,), (26) 



obtained by using Eq. (22) to calculate the state 
[H,Xjxl]\0). The first term in Eq. @ is the energy 
of the two non-interacting As, while the next two terms 
come from the Coulomb-induced LL coupling. Similar 
to \XiXi'), the state \XiY v ) = X/yJ|0) describes a non- 
interacting pair of Xi and Yv excitations. Finally, as for 
the undoped case, the last term in Eq. (|26|), 



\Bi 



= [Y?,XI]\0) = [[H,XI},XI]\0), (27) 

comes from the X-X interactions. Eq. ( |B1[) demon- 
strates that the state \Bui) is a linear combination of two 
e-h pairs with different center of mass momenta, but with 
the 2DEG in its ground state, and thus describes biex- 
citon bound, X%, and-scattering, XX, states similar to 
the undoped system£3'Lll 

The X-X and A-2DEG interactions contribute to the 
time evolution of the photo-excited 2-h state in Eq. (fTi 



7 



through |f/ , 2™*)- We further decompose the latter state 
into (a) the contribution of the non-interacting pair of 
Xj-\\pi) lh excitations, and (b) the contribution due 
to the interactions between all the different pairs of 1-h 
excitations, i.e. the X-X interactions (as in the undoped 
system) and the X interactions with the {l-/i/2DEG*} 
states (such as the four-particle Y excitations discussed 
above). 

To obtain the equation of motion of the correlated 2- 
h contribution we note that the time-evolved state 
\ip2) contributes to the third-order nonlinear response at 
second order in the applied field. By taking the time 
derivative of Eq. @ and using Eqs. @, (§§), @, 
(p5|), and (|26|), we obtain that 



P^} - p*>x\ 



l&>. 



(28) 



Recalling that \Bw), Eq. (p7|), is the interacting two- 
exciton state, we see that the first term on the rhs of 
the above equation describes thq . X pXJnteraction effects 
similar to the undoped case.ll3c30't§E2l The second term 
describes the scattering of Xi with the carriers in the {1- 
/i/2DEG*} state |^). 

Finally, we turn to the 0-h state. In Eq. ( p"8] ) we split 
this state into the contribution of the ground state |0), 
with amplitude (0\ip) = (O|-0o), and the {0-/i/2DEG*} 
contribution \tpQ nt ). The latter 2DEG* contribution is 
generated by the two-photon processes of excitation and 
de-excitation of the system accompanied by the scat- 
tering of 2DEG excitations, and is further decomposed 
into two parts. The first part, — ^ P-^* X^x) , de- 
scribes the de-excitation, after time t, of Xi from the 
{l-/i/2DEG*} state without scattering with the \-tpi) 
carriers. The latter scattering, as well as the time evo- 
lution of the 2DEG excitations created via second-order 
processes analogous to the, ones that lead to the inelastic 
Raman scattering signaled, are described by the second 
part, |-0 O > • 

The 0-h state |V>o) contributes to the third-order non- 
linear response to second order in the optical field. By 
substitutingEq. (jT|) into Eq. (JllJ) and using Eqs. @, 
(p3|) , and (|25|) , we obtain the equation of motion 



id^o) ~ H\$o) = P i* P v *i \Yi> 



P» L *Y, - P„ L *X, 



\^i)-Y, p i L * p ^ x ^ 



VS* J2 N i /2p ^ ([Xi,X$]-Sii.)\Q). (29) 



The first term in Eq. (|29|) describes the photo-excitation 
of the 2DEG via the second-order process where the exci- 
ton X^ , photo-excited with amplitude P?f , scatters with 
the 2DEG into the state IK'), and then the exciton X; is 



deexcited with amplitude P^ . The above process leaves 
the system in a 2DEG* state. It is analogous to the pho- 
toexcitation of coherent phonons in undoped semicon- 
ductors, and dominates the inelastic Raman scattering 
spectra of the 2DEG.CJ The second term on the rhs of 
Eq. (|2^) describes the scattering of Xi with the carriers 
in during its de-excitation. The rest of the terms de- 
scribe the possibility to create 2DEG excitations by pho- 
toexciting an exciton whose hole then recombines with a 
2DEG electron. 



IV. NONLINEAR POLARIZATION EQUATION 
OF MOTION 



In this section we derive the equation of motion of the 
third-order nonlinear polarization Pi{t) that determines 
the FWM and pump-probe nonlinear optical signal, and 
separate out the factorizable from the correlated contri- 
butions. When we discuss the physical meaning of the 
different terms we focus on the FWM. 

By taking the time derivative of Eq. (|9|) and using the 
definition of the operator Yj,, Eq. (|22]), we obtain that 

id t Pi{t) - n t p t {t) + Vw p i> it) = 

-fx£(t) ]T nI /2 (nx,xlu) + (1>\YM>). (30) 

The first term on the rhs of the above equation describes 
the Pauli blocking effects, which only lead to positive 
time delay FWM signal and are determined by the den- 
sity of the LLi photoexcited carriers (recall Eqs. (g) and 
(0)). The second term describes the optical signal gen- 
erated by the interactions between the recombining exci- 
ton Xi leading to the coherent emission and the photo- 
excited and 2DEG carriers. This interaction contribution 
dominates the FWM signal for negative time delays.El The 
above two source terms can be obtained by considering 
their equations of motion, which leads to an infinite hi- 
erarchy of equations of motion. Alternatively, one can 
first separate out the 2DEG* and the X-X and X-|^i) 
interaction contributions by using the decompositions of 
the photo-excited state, Eqs. ([l5|), (|l7|), and (pi), and by 
retaining contributions up to third order in the optical 
field. Using the property {0\X t ,% = {0\YiX v + (0\B H , 
(Eq. (p7|)), the expansion Eq. (|l0|), and some algebra we 
then obtain that 

= j^pt{B ii >w 2 )+j^pt'{M U ,$ )* 

i' i' 
i'j> 

+ J^pL(MYi,xtMi) + Pu (31) 

i' 

where we introduced the state 



Y\X, 



(32) 
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FIG. 3: Dominant resonant process determining the MP 
correlation contribution to the FWM signal. The first three 
panels show the Stokes Raman scattering that creates the 
MP, while the other three panels show the reverse process 
that returns the system to the ground state 



Noting that (0\Mu/) = 0, we see that the above state 
describes an excited 2DEG configuration with full valence 
band. The first term in Eq. (|3l]) describes the X—X 
interactions. Its equation of motion can be obtained by 
projecting the state (Bu>\ to Eq. (|l3|). In many cases 
it is useful to decompose the above contribution into HF 
and correlated X-X interaction contributions by using 
Eq.©: 



+ E P f ( B H'\ X j$l) + B r 



(33) 



The first term describes the ilF X-X interactions, fa- 
miliar from the undoped case.E3 The second term comes 
from the exchange process where the first optical transi- 
tion creates the {1— h/2DEG*} state and the sec- 
ond transition excites an Xj e-h pair while returning the 
conduction electrons to their ground state. The above 
process results in two e-h pairs, which scatter with each 
other while the 2DEG is at rest. Subsequently, one of the 
above pairs, i', is de-excited by the optical field, while the 
remaining pair, i, recombines and leads to the coherent 
emission. The last term in Eq. ( |33] ) describes biexciton 
and XniX. scattering correlations. Similar to the undoped 
case,E£lEZI such effects are characterized by the amplitude 
of the correlated 2-h photo-excited state, 



B«.{t) = (B u ,\,p 2 ). 



(34) 



The effects due to the propagation in time of the in- 
termediate 2DEG excitations, photoexcited via the two- 
photon process in Fig. |l](a), are described by the ampli- 
tude 



Mw{t) = (M^o). 



(35) 



Such time propagation leads to non-Mar kovian effects. 
In the case of particular interest here, the correspond- 
ing resonant contribution to the FWM signal is due to 



the nonlinear optical process shown in Fig. ||. The X 
photo-excited by the first optical transition decays into 
Y excitation. The e-h pair in this X+MP state then re- 
combines, leading to coherent emission, and leaves the 
system in a MP state. This MP propagates in time and 
then scatters with the second photo-excited X into an X 
state subsequently annihilated by the optical field. It is 
interesting to note the similarity of this process and the. 
familiar one of coherent antiStokes Raman scatterings 
that, however, involves phonons. 

The second line in Eq. ([H]) describes a shake-up of the 
2DEG during the exciton recombination that gives the 
coherent emission. In particular, the photo-excitation of 
two excitons, Xy and Xji , is followed by the recombi- 
nation of one them assisted by the shake-up of a 2DEG 
excitation. The above process leaves the system in a {1- 
/i/2DEG*} state, which is then annihilated by the optical 
field. To interpret the first term on the third line of Eq. 
( |3l| ) , we note that the HF XX interaction can be thought 
of as arising from the_scattering of the polarization with 
the coherent density.EJ Similarly, this term describes the 
scattering of the polarization with the incoherent density 
of photoexcited carriers in the {l-/i/2DEG*} state l^i). 

Finally, the last term on the rhs of Eq. ( pl| ) is the 
correlated contribution 

Pi = (V>o|o) (y#i) + J2 P 3 L *( Y * X M 



+$ |y < $i> + $i|r i $ 2 ). 



(36) 

The first two terms of Eq. (jjsTJ) describe the dephasing 
of the X and 2-X amplitudes that determine the third - 
order nonlinear polarization, while the last two terms de- 
scribe the dephasing of the incoherent contribution to 
the nonlinear polarization. Note that, by linearizing the 
above equation, we recover the correlation function P 4 - , 
Eq. (Q), that describes the dephasing of the linear po- 
larization . 

Using the above results we obtain the following equa- 
tion of motion for the third-order nonlinear polarization: 



idtPi-SliPi + ^Va'Pi' 



P 



=^£ e^ /2 (iPi\(su>-ixai])\iPi) 

i' 

"~ E (BiilXi'Xy) Pi P^.Pf* +E B «'^* 

i' 

+E^^' 

i' 

l^p£p£{o\[x f A,Y?mir 



jj'i' 



i'j' 



+E^^i|[^,^]|^)+E^' p /*^i4iv ; i)- (37) 

V ji' 

In the above equation we have separated out the source 
terms into (i) a coherent part, determined by Pauli block- 
ing effects, HF XX interactions, and the propagation 
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in time of the intermediate interacting XX and 2DEG* 
states (first three lines on the rhs), and (ii) an incoher- 
ent part, determined by the {l-/i/2DEG*} photoexcited 
state l^i) (last two lines on the rhs) and the correlated 
nonlinear contribution Pj. It is worth noting in the above 
equation that the terms proportional to the polarization 
describe a time-dependent photo-induced renormal- 
ization of the Xi energy and dephasing ( i' = i), and of 
the coupling Vw of the Xi and Xy states ( %' 7^ i). If 
the {l-/i/2DEG*} and \Ba>) excitations decay rapidly 
while the 2DEG excitations are long lived, then the non- 
Markovian LL coupling and dephasing effects are domi- 
nated by the correlation function Aiw(t). 

As demonstrated by the above equation, the Coulomb 
correlations lead to new contributions to the nonlinear 
polarization, determined by many-particle correlation 
functions. In the next section we turn to the problem 
of solving for the correlation functions on the rhs of Eq. 
([37]), and address the issue of dephasing in a strongly 
correlated electronic system. 



operators. 

100 In our case however, the ground state 
|0) may contain a strongly correlated electron gas, while, 
unlike for phonons, the MP creation operators are made 
of electrons. Thus we must use a basis of strongly corre- 
lated states that is made out of electrons. 

A basis set useful for calculating Green functions for 
tight binding and Hubbard Hamiltonians is the Lanc- 
zos basis. t£l Such correlated states can be used, to obtain 
exact solutions in the case of small systemsEJ, but also 
to approximately describe continuum resonances in large 
systems, e.g. the Fano resonances in the absorption spec- 
trum of semiconductor superlatticesEJ Each new basis 
state is obtained by acting with the Hamiltonian H on 
the previous state, and then orthogonalizing the result 
with respect to all existing basis states.E3 This procedure 
is similar to Eq. ([n]) that introduced the states \Yi), and 
led us to the parameters Sli and Vw, Eqs. (|2p| ) and (|2~l]). 
A new basis state = Zj\0) is now constructed from 
the relation 



hiya = sum 



Wiu\X, 



(38) 



V. DEPHASING AND CORRELATION 
PROCESSES 

The equations of motion for the correlation functions 
that enter in Eq. (|37]) may be obtained after introducing 
a basis suitable for describing the lh, 2h, and Oh states. 
In strongly correlated systems, an expansion in terms of 
the pair excitations of a non-interacting many-electron 
may not be convenient. In general we must introduce a 
basis of strongly correlated states that already incorpo- 
rate the ground state correlations. The choice of such a 
basis depends on the ground state and on the most im- 
portant excitations for the experimental parameters of 
interest (e.g. the filling factor, the central excitation fre- 
quency, the polarization of the optical field, etc). 

Our goal in the rest of this paper is to identify the 
dominant features in the FWM spectrum at the magne- 
toexciton energies that come from the interactions and 
time propagation of the MP collective modes. Analo- 
gous questions regarding the role of X-X interactions in 
undoped semiconductors wexe_first addressed by using 
average polarization models .BO As we show in the next 
section, a generalized average polarization model can be 
extracted from the theory developed in the previous sec- 
tions after introducing a basis of Lanczos strongly corre- 
lated states.tj This model explains the main cuialitative 
features observed in recent FWM experiments. 

We start by considering a basis for the 1-h state \4>i). 
Noting the analogy with X-phonon interactions in the 
undoped system, discussed in the introduction, we would 
like to consider a basis that consists of products of e-h 
pair and MP wavefunctions. In the undoped system, such 
states have the form e^^h^a^O) , where creates the 
phonon state, Q_ t Ji] i creates the two-particle e-h pair 
wavefunction, and |0) is annihilated by all the e and h 



where 



0,: = 



(Yi\H\Yi) 



(Yi\Yi) 

is the average energy of the four-particle excitation \ Yi 



(39) 



W iH = (Xi>\H\Yi) 



(40) 



gives the probability amplitude that Yi scatters into AV , 
and we introduced the operator 



= ^Hl-QiYi-J^Wu'^i' 



(41) 



Using Eqs. (|38|), (J39|) , and (|40|), as well as the orthogo- 
nality (Xj\Yi) — 0, one can see that the state \Zi) is or- 
thogonal to all the states \Xj), j = 0, 1, • • • , and to |Fj). 
Therefore, it is a linear combination of all the 2DEG* 
states into which \Yi) can scatter. Additional basis states 
can be constructed by applying the above orthogonaliza- 
tion procedure to the state H\Zi). We note here that the 



states \Xi),\Yd,\Zi 



do not correspond to an expan- 



sion in terms of excitations of a noninteracting many- 
electron state, since they are obtained by the action of 
the operators Xj ,Y^,zj,--- on the ground eigenstate 
|0 > of the many-body Hamiltonian H. 

By using Eq. ([lj]) and the orthogonality (Xj\Yi) = 
we obtain the useful relation 



W iH = (Y i ,\Y i )=W* % 



(42) 



Note that (Yi'\Yi) 7^ 0, and we may also have that 
(Yi'\Zi) 7^ for i' 7^ i. If this is the case we need 
to orthogonalize the independent states |^'), and then 
subtract a linear combination of the latter from \Zi) in 
Eq. ( |38| ) so that all the Z and Y states become orthogo- 
nal. However, in the electron hole symmetric limit of the 
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2DEG system, |Yj) is the same state for all i when only 
LLO and LL1 contribute (see Appendix |X|), and thus the 
latter procedure is not needed. 

Equations of motion for all correlation functions deter- 
mined by the state |V>i) can be obtained after expanding 
in a basis set of {l-ft,/2DEG*} states. Let us consider for 
example P^{t), Eq. which describes the dephasing 
of the linear polarization Pj L (t). If we choose the Lanc- 
zos basis discussed above, we obtain after multiplying 
Eq. (jH) by (Yi\ and using Eqs. ©, ©, and © the 
equation of motion 



of motion 



ir, 



(43) 



where we introduced the correlation function — 
(Zi\ipi) whose equation of motion can be obtained in a 
similar way as that of P t L . 

It is important to note that the dephasing of the op- 
tical polarization obtained as above is non-Markovian. 
Indeed, after solving Eqs. (|23|) and (fl3|) by Fourier trans- 
form we obtain that: 

[w - fi*M] Pti") + £ V w {u>)P${u>) 

= -M u )Nl<> + 2£M, (44) 
oj — Hi 

where the X energy f2j(u;) and the coupling between the 
X states Vw(lj) include frequency-dependent self-energy 
corrections due to the X-2DEG scattering, 



Q;(w) = Qi + 



-a 



o — o, 



(45) 



Additional self-energy corrections arise from Z L , dis- 
cussed in the next section. The frequency-dependence 
of the above X energies and coupling constants is a man- 
ifestation of the non-Markovian behavior of the system. 
This arises because part of the optical excitation is tem- 
porarily stored in the shake-up excitations described by 

p t- 

Using the recursive method we can also construct a 
basis for the 2-h and 0-h states, which we can then use 
to calculate the correlation functions determined by the 
states IV^) and I'tpo}- We start with the 2-h state \Bw) 
that determines the XX correlation function Ba>, and 
introduce the Lanczos state \Bw) as follows: 



H\B ii ,) = nf i ,\B w ) + \B ii/ ) , n| 



(Bw \H\Bu 
(Bw |-Bjj', 



(46) 



where Vtfj, is the average energy of the interacting 2-X 
state \Ba'). The state \B~w)> (B~w\Bw) = 0, is a linear 
combination of all the 2-X states into which \Bgr) can 
scatter. By projecting the state (Bw | to Eq. ( |28| ) and 
using Eq. (f46|), we then obtain from Eq. (34) the equation 



j'j 

+ ^P^{B iV \Yt\^)-J2^{Bw\X i j \^i)+B iV , (47) 
3 3 

where Bw(t) = (Su'lfo)- The above equation describes 
the time evolution of the "intermediate" 2-X state \Bw), 
which is created by the X-X interactions. In the case 
of undoped QW magnetoexcitons, Bw (t) corresponds to 
F(t) of Ref. 67} 

One should note here the similarity of Eq. ( [47| ) and the 
average polarization model that has been successful in 
describing the X-X correlations and biexciton effects in 
undoped semiconductors.BllSEZl This model includes the 
XX self-energy effects due to the higher Lanczos states 
\Ba'),--- via a phenomenological dephasing rate. The 
validity of such a model in the case of undoped QW mag- 
netoexcitons was analyzed in Ref. [j7"[ 

Similar to Bw, the correlation function Mw describes 
the time evolution of the "intermediate" photo-excited 
2DEG* state \Mw)- Using the Lanczos method we ob- 
tain that 



H\Mw)=n^\M v 



\Mw) , n 



M 



(M u ,\H\M n 
(Mw\Mw) 



(48) 

where is the average MP energy, and the state \Mw), 
(Mw\Mw) — 0, is a linear combination of all the states 
into which \Mw) can scatter. We then obtain, after pro- 
jecting (Mw\ on Eq. (p9|), the equation of motion 



id t M l 



Mi 



33' 



\X 3 \Y 3l )P^pt 



E P 3 L * 



(MwlYjm -Y^iMwlXjlXj,) Pf, 
3' 

-^^(Ma-I^IW 
3 

3' 3 

where Mw — (Mw\i>o) describes the dephasing of 
Aiw- In tb£_ case of the 2DEG, the single mode 



approximatiorJzj suggests that the latter dephasing can 
be treated to first approximation by introducing a phe- 
nomenological dephasing rate. 

The remaining step is the calculation of the correlated 
contribution Pi, Eq. (j36|). The equation of motion for the 
first two terms of Eq. ( |36| ) can be easily obtained from 
Eqs. ( |l2| ) and ([l3]) after using Eq. (^) and the property 

(YiXj \H = (Si + n 3 - ) (YiXj I - J2 v 33' ( Y i x j> I (50) 



■J2 w u'( x i' x 3\ + (YilYj + {Zi\Xj, 



11 



obtained by calculating the commutator [YiXj,H] using 
Eqs. (p2| ) and (|38|). The equations of motion for the last 
two terms in Eq. ( |36| ) can be obtained from Eqs. (29), 
(^||), and ( p8|) by using Eq. ( fllf ) and the basis of choice. 
We thus obtain the equation of motion 

id t Pi - n t P t -Zi=Y^ W w Pi' + Qr t , (51) 



where the correlated contribution 

Zi = (Vto|0) (Z# 1 )+^P/*(^|X i |V 2 ) 

3 

+<^ |^i> + (^ife) (52) 

has the same structure as Pi (with the difference Y — > Z;) 
and describes the dephasing of p. The factorizable con- 
tribution Qy ; describes photo-induced nonlinear correc- 
tions to the dephasing and energy of Pi , and to the scat- 
tering amplitudes Wui . The equation of motion for Zi , 
which to first order in the optical field coincides with Z(", 
has a form analogous to that of Pi . 

One should note here that the correlation function 
Pi can be decomposed further in the case of systems 
where the X-X interaction contribution to the operator 
Yi can be separated out. This is possible for example 
in undoped semiconductorsaO, where the operator Yi 
can be decomposed into a part that is independent of 
the phonon variables, which describes the X-X Coulomb 
interactions, and a part that describes the phonon cre- 
ation/annihilation processes. The former X-X contribu- 
tion comes from the last term in Eq. (]36|), and corre- 
sponds to the correlation function Z of Refs. ||,[ll] that 
mainly contributes to the six-wave-mixing spectra.u^l 
The above distinction between the interaction processes 
is possible in systems where the creation operators of 
the ground state excitations of interest (phonons, MPs, 
magnons, ■ • • ) commute with the electronic operators 
that describe the photoexcited carriers. 



VI. GENERALIZED AVERAGE 
POLARIZATION MODEL 

In this section we present an example of how the the- 
oretical framework developed so far can be used to de- 
scribe the nonlinear optical dynamics of the 2DEG in a 
high magnetic field. We consider the case where only the 
first two LLs are photo-excited, so we retain in our calcu- 
lations only the LLO and LL1 magnetoexcitons. We focus 
on filling factors close to v = 1, whem.the ground state 
2DEG populates spin-| LLO statesJ^lO and on photoex- 
citation with o+ -circularly polarized light, which excites 
spin- J, electrons.cHl The above conditions apply to the ex- 
periment of Ref. 30 that we wish to interpret. 



The electron-hole symmetry of the ideal 2D system, 
analyzed in Appendix [Al relates the correlation functions 
and interaction parameters with different LL indices that 



enter in the equations of motion. For example, in Ap- 
pendix [a| we derive the symmetry property 



(53) 



where Y is determined by Eq. (A17). The above sym- 
metry relation can be used to reduce the number of in- 
dependent variables. For example, from Eqs. ([59]), (142* 
and (pi) we obtain that = O, Ww — Wqi, 



(1 - Vi) W u = -V(1 -i/o)(1-i /i) Wqi =W= (Y\Y), (54) 
y/T^P 1 (t)=-VT^P (t) = WP(t). (55) 

where i = 0, 1, It is then convenient to make the trans- 
formation 



p -> PiVT 



(56) 



and redefine for simplicity 



Vi 



TW(1-^)(1-^) , W -> W(l - vq)(1 - vi). 

_(57) 

Using the above relations we obtain from Eq. (|2^) the 
following equations of motion for the linear polarizations: 



[fxS(t) + V 01 (l 



idtP^ 



»r )P L 



(58) 



(fio 

Vi.)P? + W(l - vi)P L ] 

idtPt = (fix - ^i)Pt 

[n£{t) + V w (l - v )P Q L - W(l - Vo )P L ] . (59) 



The above equations have the form of two coupled two- 
level systems, corresponding to the LLO and LL1 mag- 
netoexcitons. This form is due to the zero-dimensional 
confinement induced by the QW potential and the mag- , 
netic field, which leads to the novel 2DEG properties. E^lEil 
The Coulomb interactions renormalize the Rabi energy 
[i£ by a mean (local) field correction proportional to the 
polarization (analogous to the undoped systemcj), and 
by a 2DEG shake-up correction proportional to P L . 

As demonstrated by Eqs. j5q) and ([59]), the polariza- 
tion dephasing is determined, in addition to the phonon- 
induced dephasing rates Vi, by the correlation function 
P L . For weak Ti, P L dominates. In the absence of mag- 
netic field, P L describes the shake-up of FS pair exci- 
tations, and leads to a non-Mar kovian deptasing due 
to the non pcrturbative h-FS interactions BI13 For the 
experimental conditions of interest here, P L originates 
primarily from the X scattering to the continuum of 
X+MP states composed of an X and a MP with oppo- 
site momenta. An analogy can be drawn between the 
above X+MP states and the continuum^aLX+X scatter- 
ing states in undoped semiconductors .ErBEZI 

We now turn to the dephasing of the X+MP correla- 
tion function P L . This X+MP dephasing originates from 
the Coulomb-induced coupling of the Lanczos states 
\Y), \Z), \Z'),- ■ • . We note that Z L and higher corre- 
lation functions do not couple to P L , only to the ampli- 
tudes corresponding to the previous and the next Lanczos 
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states. For example, the equation of motion for Z L reads and dephasing rate 



id t Z L = 

where we defined 
(Z\H\Z) 



n z z L + w zy p l + {z'\ipi), 



ft? 



(Z\Z) 



W ZY 



(Z\Z) (Z\H\Y) 
(Y\Y) W 



(60) 



(61) 



The amplitudes of the higher Lanczos states satisfy sim- 
ilar equations of motion. After taking the Fourier trans- 
form and using the above symmetry properties and some 
algebra we obtain that 



uj - ft + ijy — Wzy^z(uj) ' 



(62) 



where 7y is the dephasing rate. The dephasing of P L is 
thus described by the self energy E^, 



Sz(w) 



1 



ft; 



%lZ 



Wz<z^z>{uy 



where Ez> is given by Eq. @ with Z' -> Z" . 

The above equation can be used to obtain a continued 
fraction expansion for the self energy. In the case of an 
TV-electron system, such an expansion terminates after N 
iterations. To obtain true dephasing for finite N, we must 
introduce the damping rates of the Lanczos states, due to 
the neglected degrees of freedom of the N — > oo system. 
The convergence of the above self-energy expansion be- 
comes more rapid with increasing damping rates.E-3 In the 
QHE literature, numerical calculations of the N-electron 
spectral functions have been extrapolated to the N — ► oo 
limit £3 

Eq. ( |63| ) can be solved analytically when the disper- 
sion in the energies and matrix elements of the higher 
Lanczos states is small as compared to the frequencies 
of interest, so that the self energy is .approximately the 
same for all the higher Lanczos states.tSl This may be the 
case for example if the momenta close to the magnetoro- 
ton minimum, q ~ l/l, dominate.cil In the case of QW 
magnetoexcitons in undoped semiconductors, the valid- 
ity of such an approximation for the X-X self energy was 
discussed in Ref. 

A microscopic determination of Sz(w) is beyond the 
scope of this paper. Here we describe the X+MP scatter- 
ing in a way similar to the average polarization model de- 
scription of the X-X scattering a In particular, we assume 
that Tiz is a sufficiently smooth function of frequency in 
the range of interest, in which case its frequency depen- 
dence can be neglected to first approximation, In the 
case of 2D magnetoexcitons in undoped semiconductors, 
a similar approximation was shown to apply for strong 
X-X interactions.EH We thus obtain the equation 



idtP 1 



(n - z 7 )p l + Pi L - Po 1 



(64) 



where the values of the renormalized Y state energy 
(Y\H\Y) 



ft 



(Y\Y) 



W Z Y~ReZ z (n) 



(65) 



7 = 1y - Wzy ImEz(ft) 



(66) 



are estimated here by-fitting to the experimental linear 
absorption spectrumEj The above approximation works 
best for sufficiently large 7. Due to the contribution to 
the Y state of finite momentum MPs and e—h pairs, we 
expect that ft > fti, where fti, Eq. ((2^), is-the energy of 
the zero momentum LL1 magnetoexciton.cH 

We now turn to the nonlinear polarization, determined 
by the equation of motion Eq.([37j). First we consider 
the X-X interaction contribution, described by the sec- 
ond line on the rhs of Eq.(|37|). For strong damping of 
the 2X states |P«/), tbcjapn Markovian X-X scattering 
effects are suppressedEaO, and we only retain the HF 
X-X interaction contribution (second term on the rhs of 
Eq.@). The X-X potentials (Bu^XjXj,), with differ- 
ent LL indices i and i' , are related to each other in the 



(63) electron-hole symmetric limit due to the property 



(1 - Ui)\Bn) = -y/(l - - Pl)\B 10 ), (67) 



where i = 0,1, that follows from Eq.(Bl). Using the 
above relation and Eq. (|56|) , we express the HF XX inter- 
action contribution in the form 



E 

a' 



yJil-vjXl-vrfiBiilXjXj,) PfPf, (P t L * - Pt) ,(68) 



where i' ^ i and the potential (Bu\XjXji) is evaluated 
in Appendix |^. 

The time-dependence of the incoherent source terms 
in the last two lines of Eq. (^) is determined by the 
{l-/i/2DEG*} state |V>i). The corresponding correlation 
functions dephase rapidly in the case of strong X+MP 
damping, unlike e.g. the correlation function Ai that de- 
scribes the time propagation of the long-lived MPs. The 
same holds for the FWM contributions due to the photo- 
induced renormalizations of the X+MP correlation func- 
tions P, Z,-- ■ , which are described by the source terms 
Qy, Qz, ■ ■ ■ in Eq. ([5l|). The latter lead to an inco- 
herent FWM contribution at the X+MP energies, which 
is broadened by the bare dephasing of P. Here we ne- 
glect such incoherent contributions to the FWM spec- 
trum. Similarly we approximate the photoexcited carrier 
density that determines the Pauli blocking contribution 
(first term on the rhs of Eq. (p^)) by the coherent exci- 
ton density, p L p L * i and neglect the incoherent density 
contribution determined by 

The correlation function Ada/ describes the time evo- 
lution of the long-lived MP intermediate states. For cr+ 
photoexcitation, the exciton operators create spin-| elec- 
trons, and thus the operators XiXj do not excite the 
spin-| 2DEG. We therefore have that 



(2DEG*\Xi\Xj 



for any excited 2DEG state. Eq. 
1. Using Eq.(p3[) and Appendix |C 



0. 



(69) 

is exact for v = 
then derive the 
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symmetry relations SI™, — Qm, -Mio — Moi, 

{l-v i ){l-v j ){M ii \X j \Y j ) = 



{l- Vl )J{\- Vj ){l-v r )(M u \X 3 \Y 3 ,) 



- -(1 - Vi )j(\- Vj ){l- Vj ,)(M ii \X j ,\Y j ) = W M , (70) 

(1 - ut)Mii = -y/(l-va)(l-vi)M i = W M M, (71) 

for any i and j ^ j'. Neglecting the incoherent X+MP 
contribution to the rhs of Eq. (|49|), determined by 
as compared to the first term, determined by the exciton 
polarizations, we obtain the equation of motion 



id t M = (n M - i lM )M 

_ r>L r>L* , pL pL* , pL pL* pL pL* 



(72) 



where the weak MP damping, which to first approxima- 
tion can be described by the dephasing rate jm ,c3o en- 
hances the non-Mar kovian dephasing effects. 

We note here that, in the absence of disorder, only 
the zero momentum MP state contributes to the non- 
linear optical signal. As already seen in the inelastic 
Raman scattering spectra, the disorder leads to the pho- 
toexcitation of a state \M) with strong contribution from 
the finite jnopentum MP's close to the magneto-roton 



energy.E^cIrcS The energy SIm is the average energy of 
the coupled MP states, and exceeds the cyclotron energy 
Q° that gives the zero momentum MP energy. 

Using the above results, and redefining for simplicity 
Wm — > Wm(1 - Vq)(1 - vi), we obtain from Eqs. fl37| ) 
and ( fill ) the following closed system of equations for the 
nonlinear polarizations: 



-2fi£(t)P L *P n L 



r>L* 



(n - *r )p - %i(i 

+W M (l - »i)M* (P L - if) - W(l - Vl )P, (73) 
id t Px = (fii - iT 1 )P 1 - V w (l - is )P Q L 

' (74) 



-2^£{t)P^Pt - 2V w (l - v )Pi L P L 



-W M (1 ~ vo)M* (P L - if) + W(l - vo)P, 

id t P=(ti-i"/)P + P 1 -P . (75) 



The second lines in Eqs. (|73j) and ([74]) describe the PSF 
effects -arid HF XX interactions similar to the undoped 
systems, while the third lines describe the correlation 
effects due to the time propagation of the intermediate 
Y and MP states. The latter correlations lead to a time- 
dependent coupling of the LLO and LL1 levels, as well as 
to non-Markovian dephasing. 

The set of four equations Eqs. @, @, @, and 
([72]), together with the linear polarization equations of 
motion Eqs. ( |5S|) , (j59|), and (p4^j, constitute our model. 
To obtain the FWM spectrum, we assume a laser excita- 
tion of the form £{t) = e^ 2 '%(i)+e* £l -%(i+Ai), where 
£ p (t) is the Gaussian envelope of the pulses emitted by 
the laser. We then solve the above equations as a func- 
tion of time t and time delay At between the two pulses, 
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FIG. 4: Linear absorption spectrum (full line). The dashed 
line shows the spectrum for P =0, in which case the dephas- 
ing is determined by the electron-phonon scattering. The 
interaction energies \/W=2.2meV and Vbi=0.5meV, the de- 
phasing rates 7 ~3meV, and the X and Y energies were cho- 
sen to reproduce the linear absorption LL peak ratio, energy 
spacing, and lineshape observed in the experiment of Ref. pOl 
for _B=8T|rMid a + circular polarization. In this case za)=0.075 
and j/i=0. E3 



keeping only the terms leading to a nonlinear signal in 
the 2/^2 — ki direction, and perform a Fourier transform 
of the nonlinear polarization to get 

P(u,At) = 0—v )P\ i (u;,At) + (l-v l )P L (u,At). (76) 

The FWM signal measured in the experiments is propor- 
tional to \P(At, uj)\ 2 and is calculated in the next section. 



VII. NUMERICAL RESULTS 

In this Section we present the results of our numer- 
ical calculations, which are based on the model of sec- 
tion VI. We start with the linear absorption spectrum, 
a(u>) oc Im[P L {uj)/£ (00)]. By fitting to the linear ab- 
sorption measurements of Ref. pQ, we can fix the inter- 



action parameters V01 and W, the energies fii and fl, and 
the dephasing rates and 7, to within ±50%. Varying 
the parameters within this fitting range yields no signif- 
icant change in the time and frequency dependence of 
the FWM or linear absorption spectrum. Fig. [| (full 
line) displays the two X peaks obtained this way. Their 
broadening is determined by (a) the X-phonon scatter- 
ing, described by dephasing rates To ~ T\ similar to the 
undoped system, and (b) the X-2DEG scattering, de- 
scribed by the correlation function P L . The important 
role of the X^X+MP scattering is clear by comparing 
to the dashed line curve of Fig. 0, obtained with P L = 0. 
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Although the X-2DEG scattering governs the lineshape 
of the LL1 peak, it plays a very small role at the LLO fre- 
quency. To interpret this behavior, we note that the main 
contribution to P L comes from {1-MP + 1-LLO-e + 1- 



LLl-h} four-particle excitations (see Eq. (A.17) and dis 



cussion m section 



Even though P L couples equally 



to both X amplitudes Pq and P±, it dominates the de- 



phasing of Pi since the above four-particle excitations 
have energy comparable to that of X±. In contrast, X 
has significantly smaller energy, and thus the broadening 
of the LLO peak is mainly determined by the X-phonon 
interactions. Note that the asymmetric lineshape of the 
LL1 resonance is due to the X+MP states and cannot be 
obtained within the dephasing time approximation. 



FIG. 5: Time delay and frequency dependence of the FWM spectrum for excitation frequency (a) at the LL1 peak, (b) shifted 
by 4meV, (c) shifted by 7meV, and (d) shifted by 9meV. The optical pulse and linear absorption spectra are displayed in the 
back panel. The pulse duration was 150fs, f2M = 17.5meV, \fWu= 2meV, and 7M=0.2meV. 



We now study the signatures of the X-2DEG correla- 
tions in the time and frequency dependence of the tran- 
sient FWM spectrum. As we discuss below, the correla- 
tion effects can be controlled experimentally by varying 
the central frequency of the optical pulse from LL1 to- 
ward LLO. This allows us to control the X amplitudes 
Pq and P\ , whose coherent superposition and interactions 
determine the FWM spectrum. Figs. || shows the effects 
of such tuning. 

Fig. ||(a) shows the FWM spectra when the optical 
pulse is centered at the LL1 peak, and the LLO peak is 
barely excited by the tail of the pulse. For such excitation 
conditions, we have that Pq <C Pi, and the photoexcited 
density of LL1 carriers far exceeds that of LLO carriers. 
As a result, the PSF and XX interaction contributions at 
the LLO energy are suppressed as compared to LL1. Dc- 

fite this however, the LLO and LL1 FWM peaks in Fig. 
a) have comparable heights. To elucidate the physical 
origin of the nonresonant LLO FWM signal, we show in 
Fig. H the contributions of PSF, XX interactions, and 
MP correlations as a function of time delay. As clearly 
seen in the above figure, for optical excitation at the LL1 
peak, the LLO FWM signal is dominated by the MP cor- 
relation contribution. At the same time, the LL1 signal 
is dominated by the PSF contribution. This is shown in 
Fig. @. The origin of the strong MP correlation contribu- 
tion to the LLO FWM signal can be seen by comparing 
the latter for different values of the MP energy Um , while 
keeping the rest of the parameters constant. As demon- 
strated by Fig. |^, the LLO signal, dominated by the 
MP-mediated LL coupling due to the process of Fig. [|, 
is resonantly enhanced as CIm approaches the Xq — > X\ 
excitation energy (~ 18meV here). 

One should note here that the XX interactions also 
couple the two LL's. However, the corresponding LLO 
signal is weaker due to the absence of a resonance, sim- 
ilar to the MP correlation signal for nonresonant Ojvf, 
and cannot fully account for the strong LLO signal ob- 
served in the experiment of Ref. To see this, note 
that, in the undoped system, where only the XX inter- 
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FIG. 6: MP correlation (dashed line), PSF (dotted line), and 
XX interaction (dashed-dotted line) contributions to the full 
FWM signal (full line), calculated at the LLO peak frequency 
for photoexcitation as in in Fig. ^(a). 



actions contribute, the FWM signal at the LLO energy 
is negligibleEa. More importantly, in the experiment of 
Rcf. |30|, the LLO peak in the doped system was sup- 
pressed as compared to the LL1 peak as the density of 
photoexcited carriers approached that of the 2DEG. In 
this case the MP correlations of the cold 2DEG diminish, 
and the doped and undoped QW FWM signals start to 
look similar .El 

We now turn to the temporal profile of the FWM sig- 
nal. Fig. ^, which plots the normalized time-dependent 
LLO and LL1 signals, demonstrates the difference in the 
dynamics between the PSF and MP correlation effects 
that dominate the LL1 and LLO signal respectively. As 
already known from undoped semiconductors, the Pauli 
blocking effects cannot lead to a FWM signal for neg- 
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Time delay (fs) 

FIG. 7: PSF (dotted line) contribution to the FWM signal 
(full line) at the LL1 peak frequency for photoexcitation as in 
Fig. |B](a). The XX interaction (dashed-dotted line) and MP 
correlation (dashed line) contributions are negligible here. 




Time delay(fs) 

FIG. 8: FWM signal at the LLO peak frequency for different 
values of the average MP energy: f2M=14mev (dotted line), 
17meV (full line), 18meV (dashed line), and 20meV (dashed- 
dotted line). Photoexcitation conditions as in Fig. H(a). 



ative time delays. The XX interactions lead to such 
a FWM signal, with rise time determined by the de- 
phasing of the interacting XX staie \B)B Within the 
time dependent HF approximation^ the latter rise-.time 
is ~ (2To + 2ri) _1 in the case of interest hereJij For 
r ~ Ti, this rise time is about one half of the XX FWM 
decay time, ~ (2ri) _1 , or the PSF FWM decay time, 
~ (2To) _1 at the LLO energy. Fig. || however shows an 
almost symmetric temporal profile of the LLO FWM sig- 



nal, unlike for the LL1 signal. The latter is dominated by 
the PSF contribution, and is thus suppressed for negative 
time delays, while the LLO signal is dominated by the MP 
correlations, and is enhanced for negative time delays, 
similar to the experiment of Ref. The origin of this 
time dependence can be seen from the equation of mo- 
tion of the MP correlation FWM source term in Eq. (|7^) . 
After retaining only resonant terms, one can see that the 
rise of this signal is governed by the time dependence of 
the product P^Pq, while the decay is determined by the 
time dependence of P±, A4, and by quantum interference 
effects. As discussed above, due to the X-2DEG scatter- 
ing, Pi dephases much more strongly than Pq. Thus the 
time dependence of PiPq is similar to that of P^, which 
results in an almost symmetric FWM temporal profile at 
the LLO frequency. Furthermore, the quantum interfer- 
ence and beating effects enhance the decay of this signal 
for positive time delays. 

The relative magntitude of the MP correlation versus 
the PSF/XX mean field FWM signal can be controlled 
experimentally by changing the central frequency of the 
optical pulse. Fig. ||(b) shows the time-dependent FWM 
spectrum for excitation conditions such that \Pq\ ~ \P\\- 
The LLO signal now dominates, and retains a temporal 
profile similar to Fig. |^(a). Note that, due to the in- 
creased pulse overlap with LLO, the PSF source terms 
in Eqs. (73) and ( |74| ) now have comparable magni- 
tude, while the XX interaction and MP correlation source 
terms are also enhanced. However, the strong dephas- 
ing of Pi discussed above suppresses the LL1 FWM sig- 
nal. This effect is magnified in the nonlinear spectra as 
compared to the linear absorption. Importantly, due to 
the resonant enhancement of Fig. Bl the magnitude of 
the MP correlation FWM contribution is enhanced more 
strongly by the increased pulse-LLO overlap as compared 
to the mean field FWM signal. Finally, weak oscillations 
as a function of time delay, with a period equal to the 
spacing of the two LL peaks, start to appear. 

As can be seen in Figs. ||(c) and[|(d), the oscillations of 
the LLO FWM signal, as function of time delay, become 
more pronounced as the optical excitation frequency is 
shifted from LL1 toward LLO. It is important to note 
that the LL1 signal is almost completely suppressed, es- 
pecially for the excitation frequency of Fig. |^(d), and 
therefore there are no significant oscillations in the real 
time t, i.e. the time related to the frequency uj via Fourier 
transform. Thus the oscillations observed in Fig. ||(d), as 
well as in the experimental datAjOiRef. ||(], have a strong 
quantum kinetic contributionJjOLj The physical origin of 
such an effect can be seen by plotting in Fig. [Io| the PSF 
and MP correlation contributions to the LLO signal as 
a function of time delay for photoexcitation as in Fig. 
||(d). PSF leads to negligible oscillations, while the MP 
correlation leads to strong oscillations. To see the origin 
of the latter, we also plot in Fig. [l(] the MP correlation 
signal obtained after neglecting the LLO coherent density 
source term PqPq* in the equation of motion Eq. ( |72] ) 
of M(t); the oscillations diminish in the latter case. To 
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FIG. 9: Temporal profile of the FWM spectrum at the LLO 
peak frequency (full line) and the LL1 peak frequency (dashed 
line). The two signals have been normalized for clarity. 
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FIG. 10: PSF (dashed-dotted line) and MP correlation (full 
line) contributions to the FWM signal at the LLO frequency, 
normalized to unity for clarity, for photoexcitation as in Fig. 
^|(d). We also plot the MP correlation signal (dotted line) 
without the LLO coherent density source term of M(t) in Eq. 
(|r2|). Note the almost complete absence of oscillations in the 
latter and in the PSF contribution. 



interpret all this, we note that, for the excitation con- 
ditions of Fig. ||(d), we have that P L > P X L , and the 
density of LLO carriers far exceeds that of LL1 carriers. 
In fact, here the PSF contribution exceeds the MP corre- 
lation contribution. Most importantly, the LLO coherent 
density source term of M. is now larger than the source 
term P^Pfr* that gives the resonant MP contribution. 



Even though PqPq* gives a nonresonant contribution to 
Ai (t) , as Pq exceeds P X L this contribution becomes com- 
parable in magnitude to the resonant contribution due 
to P^Pq* . The beating between the above two resonant 
and nonresonant processes gives rise to the strong oscilla- 
tions. By shifting the excitation frequency further toward 
LLO, eventually the PSF contribution dominates, and the 
FWM dephasing is determined bjj— the electron-phonon 
and intra-LL dephasing processes. Ej 



VIII. CONCLUSIONS 

In summary, we presented a theory that provides a uni- 
fied description of the ultrafast nonlinear optical response 
of a large class of semiconductor systems with a strongly 
correlated many-electron ground state. Our main result, 
Eq. ([}?]) , gives the equation of motion for the third-order 
nonlinear polarization measured in transient wave mix- 
ing and pump-probe experiments, and allows us to study 
the role of the correlations and the interplay between co- 
herent and incoherent effects. Our expansion in terms 
of the optical field is valid for sufficiently short pulses 
and/or weak excitation conditions, where the correlations 
are most pronounced. Our theoretical framework allows 
us to describe the role of the long-lived collective excita- 
tions of a strongly correlated cold electron gas, which is 
present prior to the optical excitation. 

Our theory was applied to the case of the 2DEG in a 
strong magnetic field. Our numerical solution for pho- 
toexcitation close to the LL1 energy with o + — a + cir- 
cularly polarized light suggests new experimental signa- 
tures of collective and correlation effects. In this case 
the relevant 2DEG collective excitations are the long- 
lived inter-LL magnetoplasmons, which dress the pho- 
toexcited magnetoexcitons and lead to polaronic-like ef- 
fects and strong non-Markovian dephasing. We showed 
that such effects dominate the time delay and frequency 
dependence of the transient FWM spectrum. FWM spec- 
troscopy using femtosecond optical pulses provides both 
the time and the frequency resolution necessary to access 
this new regime of 2DEG physics. Our theory allowed us 
to study in a systematic way the experimental signatures 
of the 2DEG quantum dynamics. We predicted, in par- 
ticular, a resonant enhancement of the lowest LL FWM 
signal, a strong dephasing of the next LL magnetoex- 
citon, a symmetric FWM temporal profile, and strong 
oscillations as function of time delay with a strong quan- 
tum kinetic contribution__Such predicitions agree with 
recent experimental dataH3 

The above correlation-induced dynamics can be con- 
trolled by tuning the central frequency of the optical ex- 
citation between the two lowest LLs, which changes the 
coherent admixture of the two MP-dressed magnetoex- 
citons, or via coherent control experiments using phase 
locked optical pulsesB Such experiments, as well as a~ 
circularly polarized optical pulses, provide new ways to 
access the very early dynamics of the strongly corre- 
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lated 2DEG, during time scales shorter than the duration 
of the interactions. Such temporal and spectral resolu- 
tion opens up new ways to observe fractional QHE non- 
instantaneous correlations, as well as magnon, exciton- 
magnetoroton, charged exciton, and skyrmion effects. 
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APPENDIX A 

In this Appendix we derive some useful expressions 
for the operators Yl, Eq. (|22|), in the case of the ideal 



2D system displaying electron-hole symmetry. To de- 
scribe the magnetic field effects, we choose to work in 
the Landau gauge A = (0, Bx, 0). The eigenstates of the 
kinetic energy operator are then characterized by the y- 
component of the momentum, k, and the LL index, n. 
The electro n, i p^ , and hole, ip a , eigenstates in this gauge 
are given byc30 

V>a(r) = ~Jf^n{x - x k ) , xj; a {r) = Ca(r), (Al) 

where a = (k,n,cr), ~a = (— k,n, cr), and the spin-cr 
wavefunction is kept implicit. In the above equation, 
fy n is the eigenfunction of the ID harmonic oscillator 
with frequency equal to the cyclotron frequency, x k = kl 2 
is the x coordinate of the cyclotron orbit center, / = 
(Hc/eB) 1 / 2 is the magnetic length (Larmor radius), and 
L is the system size.B^EZl 

The operator Yi is determined by the commutator 
[Xi, H int ] , where the Hamiltonian H rnt = V ee + V h h + V eh 
describes the Coulomb interactions: 



H, 



dvdr' v(r - r') ^(r)^(r) - ^(r)ip(r) ip* (r')i/>(r') - ^ (r')<if)(r') 

I 



(A2) 



where (r) is the electron creation operator, %j> t (r) is the basis wc transform the Hamiltonian Eq. ( |A2| ) into the 
hole creation operator, and v(r) is the Coulomb potential. familiar form 
By expanding the above creation operators in the Landau 
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(A3) 



where, in the ideal 2D system, the Coulomb interaction where v q = 2ire 2 /q is the Coulomb potential, and 
matrix elements v l d ia2 Q3Q4 (with i,j = e,h) are given by 



aia2,c*3CK4 



/ ^Vl Q2 (q)*<U(-q), (A4) 



K ia M) = / dvr ai {vY*^ a2 (v) , i£ aa (q) = J ^^We^-^W. (A5) 



r 



Following Ref. ^6] we obtain that where 

K,«M)=Vn 1 nM)h,kM^ (A6) fk lk M=e tq * {kl+k2)l2/2 S kl M+ qy (A7) 
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and, for m > n, we have that 
(-q v + iq x )l 



^mn(q) = —7 

m! 



V2 



m-n / 2/2 

m-n I V/ 1 



V: 



(A8) 

where L™ n is the generalized Laguerre polynomial. 
<Pmn(q) for m < n can be obtained by using the property 



yw(q) = 92* TO (-q). 



(A9) 



Using Eq. QAlh we obtain from Eq. (A5) 



*&a,(q) = (q)- 



(A10) 



The following symmetry relations can be shown by using 
the above relations: 



v- = v 31 v eh = v ee v hh = v ee 

^CKl Q2 ,Q304 ^0.30.4,0110.2 ' 0102,-04 — 03 u QiQ2,tt3Q4 ' 04 — 03, Oi O2 ^0304,-02-^1 " 



(All) 



The com mutator \h- a e ai H j nt ] can be calculated from 
Eq. (A3). Using Eq. (All) and some algebra we obtain 



that 



[h- a e a , H int ] = - v 



Oi Q2Q' 



\v ee ,(e ] t 

u ai 02 ,oa' \ oi L 



h J h \h p , 

l(/ —a.2 a i 1 Q a 



(a <-> a') . (A12) 



After summing over k, and recalling the definition Eq. and 
([|) of the X operators and the definition of N na , the 

lhs of the above equation becomes the commutator 

1/2 ~ 

NJa [X na ,H mt ]. Using the properties 



/,fkk 2 (q)/fcifc(-g) = 4i* 



*I u(g)<^„n 3 (q) </> ni „(- q) = 5 nu „ 2 \ dqv(q) \(j> nni (q)| 2 , 

(A14) 



(A13) we obtain after using Eq. (|4|) and some algebra 
I 



[X n(T , Hi n t] — — K?n'o-(l — v n'a)^-n'a 



n) £ [C ia2 ,M'n'^™*"''"( Tif4,i ') ( A15 ) 



r 



where 



V° , = 



1 



7^2 «9 |0n»'(g)| 2 - 



•\/(l - l-w)(l - V n i a ) 

(A16) 

We now restrict to the first two LL's, which dominate the 
optical spectra for the excitation conditions of interest. 
Recalling Eq. ( |22| ) we see that the operator Y n is deter- 
mined by the last term of Eq. (A15). The only nonzero 



contribution to this term comes from n' ^ n, and there- 
fore n'=l if n=0, or n'=0 if n=l. As a result, the rhs of 



Eq. (A15) changes sign between n=0 and n=l, and we 
obtain Eq. (|53|). The explicit expression for the opera- 
tor Y a = yT — v\ a Y\ a can then be obtained straightfor- 
wardly by subtracting the X contributions defined in Eq. 
(p2|) from the operator 



J 



pp' kk' mm' a 



(^mfT'V" 1 '* 7 ' ^-p'mV'^-P^'J (^pmp'm , ,klk / 0^ l —kla^k'0(X v pmp / m / , kOk'l^—kOo-^k'laJ ■ (A17) 
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The subtracted X contributions describe corrections to 
the X energies and Coulomb-induced LL coupling due 
to the 2DEG. As discussed in section [vj, for photoex- 
citation with tr + circularly polarized light, we have that 
a =|. For filling factors close to v = 1, the spin - J. states 
are empty. We can then decompose Eq. ( A17 ) into the 



a 1 =1 term, which describes the X-X interactions, and 



the a' 
tions. 



-1 term, which mainly describes X-MP interac- 



APPENDIX B 



In this Appendix we evaluate the HF XX potentials 
(Bu\XjXji) in the ideal 2D system. We consider tr + pho- 
toexcitation and filling factors close to v = 1 so that Eq. 
(|69| ) applies. Recalling the definition Eq. (^) we obtain 
from Eq. ( A17 ) after using the property [Xi,Xj] = 
that 



{B\ a ',ncj\ — 7= = 5^ \ V pnp'n',klk'oh 

y/N na N la , pp , kk , n , L 



-pna&p' n 1 a^ — kla' &k'0a' '^pn / p'n,klk'0^'~ pn' a^-p' na^— klcr' ^fc'Oo 



v pnp ! ri ,k$k' \ h-pna&p' ii' ah-kOa' &k' la' + V p n / p'n,kOk f 1 h-pri a&p' na h-kOa' &k' 



(Bl) 



The only nonzero contribution to the above equation 
comes from n' ^ n. Noting the LL indices, we see that, 
for the conditions considered here, we have that 



(Bu\XjXj) 



0. 



(B2) 



Substituting the definition of Yi, Eq. (|22|), into Eq. (|27|), 
restricting to the first two LLs, and denoting i' i, we 
obtain that 



APPENDIX C 

In this Appendix we derive some useful relations for 
the overlap (2DEG*|M^), where \2DEG*) is any 2DEG 
excited state, for filling factors close to v = 1 and for tx + 
polarized light. Using Eqs. (p3) and (53) we obtain that 



(BulX.Xa) = (X^lHXllXo) - (X^HX^X.Xo) 
-SliiXiX^Xo) + VwiXtXolXtXo) - (YHXilXtXo^Bd) 

We have that (XiXi\XiX ) = {XiXi\X v X v ) = due 
to the orthogonality of the valence hole states, while 
{Y^X^XxXq) = due to Eq. ©. Using the above, 
Eq. (||) for the commutator [H, xjl Eq. (|l|) for the 
states H\Xq) and (Xi\H, and Eq. (B2), we obtain after 
some algebra that 



iv 1/2 |M 0l ) = iv 1/2 y A;|o) 

= -N^YXjlO) = -Nl /2 \M U ). 
From Eq. (p2) we obtain after using Eq. (p2|) that 



(CI) 



\M W ) = Xi\Y v ) -{H + Qi- Cl^XilXi,) 

+ J2 VijXjlXi,) - £ v ji'Xi\Xj). (C2) 



(BalXiXo) = 2V ii >(X 1 X \X 1 X ) 
-VoxiXiXilXoXo) - VioiXiXilXtXt). 



Using the relations (XoXi\XiXq) = 1 and 

(XiXilXiXi) = 2(1 
obtained from Eq. (^|), we finally obtain that 



.1) 



N i (B ii \X 1 X ) = 2V ii , , i'^i. 
The above relation recovers the results of Ref. 43 



(B4) 



(B5) 



(B6) 



The state Xi\Yii) describes a 2DEG excitation, created 
via the process shown in the first three panels of Fig. (|^). 
Using Eq. @ and the property (2DEG*\H\0) = we 
obtain that 



{2DEG*\M iV ) = (2DEG*\X i \Y i , 
Using Eq. (p3|) we then obtain that 



(C3) 



7V 1/2 (2DEG*|M i0 ) = No /2 (2DEG*|Xi|Y ) 
-iV 1 1/2 (2DEG*|X i |Yi> = -N} /2 (2DEG*|Mii). (C4) 
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